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ABSTRACT 

We survey the properties of all orbit families in the rotating frame of a family of realistic triaxial potentials 
with central supermassive black holes (SMBHs). In such galaxies, most regular box orbits (vital for maintaining 
triaxiality) are associated with resonances which occupy two-dimensional surfaces in configuration space. For 
slow figure rotation all orbit families are largely stable. At intermediate pattern speeds a significant fraction 
of the resonant box orbits as well as inner long-axis tubes are destabilized by the "envelope doubling" that 
arises from the Coriolis forces and are driven into the destabilizing center. Thus, for pattern rotation periods 
2 x 10 8 yr < T p < 5 x 10 9 yr, the two orbit families that are most important for maintaining triaxiality are highly 
chaotic. As pattern speed increases there is also a sharp decrease in the overall fraction of prograde short-axis 
tubes and a corresponding increase in the retrograde variety. At the highest pattern speeds (close to that of 
triaxial bars), box-like orbits undergo a sudden transition to a new family of stable retrograde loop-like orbits, 
which resemble orbits in three-dimensional bars, and circulate about the short axis. Our analysis implies that 
triaxial systems (with central cusps and SMBHs) can either have high pattern speeds like fast bars or low patten 
speeds like triaxial elliptical galaxies or dark matter halos found in A/-body simulations. Intermediate pattern 
speeds produce a high level of stochasticity in both the box and inner long-axis tube orbit families implying 
that stable triaxial systems are unlikely to have such pattern speeds. 



1. INTRODUCTION 

It is widely accepted that since both elliptical galaxies and 
dark matter halos form via hierarchical mergers, they should 
be triaxial. If a significant amount of gas is present in the 
progenitor galaxies, its dissipative condensation generally re- 
sults in systems that are more oblate than those produced 
in purely collisionless collapse ([Du binski 1994, Kaza ntzidls| 
|et al. ||2004) |Debattista et al.||2008| l However, recent stud- 
ies have shown that two orbit families that characterize tri- 
axial systems (boxes and long-axis tubes) persist in signifi- 
cant numbers even when a triaxial system looks al most oblate 
( |van den Bosch & de Zeeuw|201 0H Valluri et al-pOlO) . Since 
merging systems generally have angular momentum (either 
internal angular momentum of the individual progenitors or 
the angular momentum of their relative orbit or both), merger 
remnants generally have angular momentum. If the remnants 
are triaxial, they can exhibit figure rotation: a property that is 
independent of the streaming motions of individual particles. 

Although it is logical to assume that triaxial elliptical galax- 
ies can have figure rotation, it is currently difficult (if not im- 
possible) to observationally distinguish between figure rota- 
tion (tumbling) and orbital streaming. Although figure rota- 
tion was first proposed to explain "anom alous dust lanes" in 
triaxial elliptical galaxies (van Albada et al. 1982), there have 
been only a few observational attempts to measure the pattern 
speeds of either elliptical galaxies o r dark matter halos ( |Bu-| 
|reau et al.|[T999] |Jeong et al ; 1[2007] l. Cosmological AMxxry 
simulations (without dissipation) predict that dark matter ha- 
los are significantly triaxial and that a majori ty of dark matter 
halos (—90%) have figure rotation ( jBailin & Steinmetz|2004 



Bryan & Cress||2007| l with slow pattern speeds following a 
lognormal distribution centered on VL p = 0.148/zkm s" 1 kpc -1 
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Bailin & Steinmetz] ( |2004[ ) found that the pattern speed of 
the figure rotation is correlated with the cosmological halo 



spin parameter A (Peebles 1969), but is independent of halo 
mass. Figure rotation of triaxial dark matter halos has been 
suggested as a mechanism for driving spiral structure, warps 
and bars in spiral and elliptical galaxies with low mass disks 
( Bureau et al.|1999||Dubinski & Chakrabarty|2009] l. The only 
family of triaxial stellar systems with strong figure rotation 
is the bars of spiral and lenticular galaxies. Rapidly rotating 
(triaxial) bars can be almost completely regular (Pfen niger & 
Friedli ]| 1 99 1 ]> and have a compl ex orbital s tructure (|S kokos 
et al. 2002a b i Patsis et al. 2002; Harsoula & Kalapotharakos 
2009). 



Emsellem et al. (2007) introduced a luminosity weighted 
measure of the specific line-of-sight angular momentum, A^, 
to quantify the angular momentum content of E/S0 galaxies 
observed in the SAURON sample. They found that this pa- 
rameter allows E/S0 galaxies to be classified into two sub- 
groups - the "fast rotators" with A« > 0.1 and "slow rotators" 
with \r < 0.1. The slow rotators are likely to be triaxial el- 
lipticals since they frequently display isophotal and kinematic 
twists, including kinematically decoupled cores. T hey tend to 
be more luminous and have shallo w central cusps (Emsellem 
et al.|2007||Cappellari et al.|2007) . In contrast, "fast rotators" 



are more numerous and close to axisymmetric, but may retain 
a significant population o f triaxial o rbit fami lies such as box 
and long-axis tube orbits ( jvan den Bosch & de Zeeuw] 2010). 
These authors find that even a small fraction of box-like orbits 
can significantly alter the dynamical estimates of supermas- 
sive black hole (SMBH) masses. Recent A/-body simulations 
of mergers of gas rich disk galaxies that include star forma- 
tion and dissipation show that both fast and slow rotators can 
be formed in mergers between disk galaxies. The latter class 
primarily form in maj or mergers with s maller gas fractions 
and are more triaxial ( |Jesseit et al.||2009| l. (Recently, |Bois| 
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|et al.| ( |20T0) l have shown that due to numerical resolution ef- 
fects in A^-body simulations, the fraction of merger remnants 
that are triaxial slow rotators may have been underestimated). 
Thus despite the absence of concrete observational measure- 
ments, some fraction of triaxial elliptical galaxies and dark 
matter halos could have figure rotation. 

There has been little work on the effects of figure rotation 
on triaxial galaxies in the last two decades - none o n realis- 
tic den sity profiles with central cusps and SMBHs. Binney 
( | 198 1 [ > studied the effect of figure rotation on the behavior of 
the closed periodic orbits in the equatorial plane perpendicular 
to the axis of rotation (z) of flattened galaxy models including 
those with weak bars. He showed that closed retrograde or- 
bits in the x-y plane lying in an annular region (now called the 
"Binney instability strip") become unstable to perturbations 
perpendicular to the plane. He also showed that this instabil- 
ity strip moved inward as the speed of rotation of the figure 

in creased. 

Heisler et al. ( 1982) studied the stability of closed orbits in 
a triaxial model with a central density core subjected to both 
slow (pattern rotation period T p ~ 10 9 yr) and fast rotations 
(T p ~ 10 8 yr). They showed that the 1:1 periodic orbits that 
circulate about the long (x) axis of the model are stable to 
figure rotation and are tipped about the y-axis by the Cori- 
olis forces in a direction that depends on their sense of mo- 
tion. Two such orbits exist: one rotates clockwise about the 
x-axis and the other circulates anti-clockwise. However, both 
"tipped" orbits circulated about the short axis in the s ame di- 
rection. (These orbits were termed "anomalous" by van Al- 



bada et al. ( 1982| l.) The long-axis tube family that is "par- 



ented" by the anomalous orbits was therefore also expected 
to be stable and also "tipped" at an angle. They noted that 
for orbits with very large energies (i.e., extending to large 
radii) such orbits could be tipped by ~ 90° into the equato- 
rial plane. Our more general analysis in models with realistic 
density profiles shows that while outer long-axis tubes indeed 
behave as predicted, the inner long-axis tubes with the small 
peri cent ric radii are easily destabilized by figure rotation (Sec- 
ti on^. 

Heisler et al. ( 1982) also showed that the normal retrograde 
orbits that he in the equatorial plane were stable (except in 
the Binney instability strip, where they were unstable to ver- 
tical perturbations). The sequence of closed, stable period or- 
bits identified by them were found to exist at both slow and 
fast pattern speeds and constituted one composite sequence 
which was stable over the entire energy range. They specu- 
lated that this implied that triaxial galaxies could have both 
fast a nd slow pattern speeds. 

|de Zeeuw & Merrltt| ( |1983| l complemented this work, with 
a study of orbits in the principal planes of a rotating triax- 
ial galaxy with a central core. They found three prograde 
sequences of stable orbits in addition to the retrograde se- 
quences. Inside the core of the galaxy, they found that the 
x-axial orbit was stable and generated a family of box orbits 
that underwent "envelope doubling" as they looped around the 
center due to the Coriolis forces. 

The first attempt to construct a self-consistent triaxial 
galaxy with slow figure ro tation was presented in the pioneer- 
ing work of Schwarzschild ( 1982 1. Based on previous studies 
of periodic orbits, Schwarzschild restricted his rotating triax- 
ial models to a moderate pattern speed (T p ~ 1.2x 10 9 yr) for 
which the main resonances (corotation, outer inner Lindblad, 
and the Binney instability strip) lay outside the model. He 



showed that self-consistent triaxial models with figure rota- 
tion could be constructed such that they had either no stream- 
ing motions or maximal streaming. 

The main conclusion from the early work on figure rotation 
in triaxial galaxies with cores was that except for the Bin- 
ney instability associated with retrograde orbits confined to 
the equatorial plane, several stable sequences of orbits that 
parented the major orbit families existed over a range of pat- 
tern speeds and energies, implying that triaxial galaxies with 
both fast (rotation period T p ~ 10 6 yr) and moderate rotation 
speeds (T p ~ 10 9 yr) could exist. Our findings here will show 
that the presence of a central cusp or SMBH alters that con- 
clusion; 

|Gerhard & Binney] ( |1985| l first proposed that the box or- 
bits which form the back bone of triaxial elliptical galaxies 
would become chaotic due to scattering by the divergent cen- 
tral force arising from a central black hole. The presence of a 
significant fraction of chaotic orbits results in chaotic mixing 
which can cause secular r elaxation of orbits in phase space 
(Kandrup & Mahon 1994) resulting in a change in the shape 
of the galaxy from triaxial to axisymme tric on timescales of 
order ~ 30-50 dynam i cal times (Merritt & Valluri|1996[|Mer-| 



ritt & Quinlan|1998| l. [Gerhard & Binney| ( |1985| l had also ar- 
gued, however, that if the triaxial potential had figure rota- 
tion, then box orbits (which are crucial to maintaining triaxi- 
ality) would be less affected by the central force in the rotating 



frame due to "envelope doubling" (de Zeeuw & Merritt 1983 ). 

Several early studies indicated that triaxial galaxies with 
figure rotation could often hav e significant fractions of chaotic 
orbits The studies by (|Udry & Pfenni ger||1988[ |Martinet & 



|Udry|1990| |Udiy||1991[|Tsuchiya et al.|1993| ) were however 
restricted to limited numbers of orbits in models with cores 
and did not pursue the primary cause of this chaos. 

It is now believed that all elliptical gala xies have either shal- 
low or steep central density cusps (Gebhardt et al. 1996 |Lauer 



eTaT1[2007ll and cent ral SMBHs (Ferrarese & Merrittj|2UUU| 



Gebha rdt et al.|[2000] >. In realistic elliptical galaxy models 



with cusps and black holes, a large fractio n of phase space 
is occ upied by resonant and chaotic or bits (Miralda-Escude 
| & Schwarzschild|1989HSchwarzschild|1993l|Valluri & Mer- 



ritt|1998||Merritt & Valluri|1999||Poon & Merritt 2004}. Al- 
though the fraction of chaotic orbits increases with an in crease 
in the strength of the density cu sp ( Sch warzschild|199"3| |Mer-| 
|ritt & Fridman|1996[|Meiritt| 19971) or the mass of the central 
black hole (Merritt & Valluri 1996), it is still possible to con- 
struct triaxial models that do not evolve rapidly due to chaotic 
mixing so long as a significant fract ion of orbits (~ 50%) are 
regular ( [Poon & Merritt|2002l[2004| . 

In a precursor to the present paper, |Valluri| ( |1999| ) first 
showed that box orbits in triaxial galaxies with cusps and 
black holes are destabilized by moderate amounts of figure 
rotation because the envelope doubling acts to further desta- 
bilize the resonant box orbits, rather than stabilize chaotic or- 



bits as predicted by Gerhard & Binney ( 1985 ). This paper 
further explores the cause of the destabilization of box-like 
orbits and investigates the effects of figure rotation on all ma- 
jor orbit families. 

Our objective in this paper is to study the behavior of all 
the major and minor orbit families, especially those that form 
the backbone of realistic rotating triaxial galaxy models. The 
goal is to identify the range of rotation speeds for which such 
orbits will remain stable. 

The paper is organized as follows. In Section[2j we describe 



Triaxial Galaxies with Figure Rotation 



3 



the numerical model, the selection of orbital initial conditions, 



and the Laskar mapping method (Laska r[l990| l. In Section[3j 
we describe the results of our analysis of the effects of figure 
rotation on major orbit families, as a function of pattern speed, 
orbital energy, and the shape of the triaxial model, in models 
with and without central SMBHs. We summarize our results 
and discuss the implications of our findings in Section|4] 

2. MODELS AND NUMERICAL METHODS 

2.1. Density model 

We studied tri axial gener alization s of spherical model s first 
presented by Dehnen|( |l993| l andlTremaine et al. ( 1994J. The 
model, which we will henceforth refer to as the "Dehnen- 
model", has a density law that is a good fit to the observed 
luminosity profiles of ellipticals and the bulges of spirals, and 
is given by 



p(m) ■■ 



(3-7)M 
4irabc 



m 7 (1 +m) 



-(4-7) 



0< 7 <3, 



where 



m = + 



b 2 



a>b>c>0, 



(1) 



(2) 



and M is the total mass of the model. The mass distribu- 
tion is stratified on concentric ellipsoids where a,b, and c 
are the scale lengths of the semimajor, semi-intermediate, and 
semiminor axes of the model, which are aligned with Carte- 
sian coordinates x,y, and z, respectively. The parameter 7 
which determines the logarithmic slope of the central density 
cusp ranges observationally from 7 = 0.5-1 in luminous galax- 
ies with "shallow cusps" (sometimes called a "core ") to 7 = 2 
in the lower luminosity ga laxies with "steep cusps" ( Gebhardt 
|et al.|199"6}|Lauer et al.|2007] >. At large radii the density pro 
file of the model always falls as m~ A . The model has a fi- 
nite density core for 7 = and an infinite central density for 
7 > 0. The potential ($(x)) and forces in the stationary tri- 
axial Dehnen model in ell ipsoidal coordinates are taken from 
|Merritt & Fridman| ( [T996] l. 

In the rotating frame, the energy of an orbit is not an integral 
of motion but the Jacobi integral (Ef) is a conserved quantity: 



i|i| 2 + <i>-I|^xx| 2 , 



(3) 



where x and x are three-dimensional spatial and velocity vec- 
tors, respectively. 
In our models, figure rotation is about the short (z) axis, 

hence fl p = fl p e z and equations of motion in the rotating frame 
(BT08, Section 3.3.2) become 

x = -V$ - 2(Q P X x) - Q p X (Q p x x) (4) 

= -v$-2(n p xx)+\n p \ 2 x. (5) 

In Cartesian coordinates the equations of motion are then 
given by 



<9$ 2 
.v = — - — 2il p y+ilpX, 



dx 

d>S> , 

y=- —+2n p x+np, 



z=- 



<9$ 

3? 



(6) 
(7) 
(8) 



where —2Q, p y and 2Q, p x are Coriolis force terms and tt 2 x and 
Q 2 y are centrifugal force termsr 

Most of the models used in this study are close to max- 
imally triaxial with triaxiality parameter T = (a 2 - b 2 ) / (a 2 - 
c 2 ) = 0.58 and minor to major axis ratio c/a = 0.5. For a lim- 
ited number of models, we also explored the effect of chang- 
ing the shape of the triaxial figure with c/a = 0.5, 0.8, and 0.7 
with triaxiality parameters T = 0.1,0.9, and 0.3, respectively. 
The shapes therefore range from nearly oblate, through oblate 
triaxial to nearly prolate. 

Following the standard practice, we adopt a set of units 
where the total galactic mass M, the semimajor axis scale 
length a, and the gravitational constant G are set to unity. 
When a central black hole is added to the model, its mass Mbh 
is expressed as a fraction of the total galaxy mass M. In this 
paper we restrict ourselves to studying models with either no 
central point mass (Mm = 0) or models with = 0.001 (i.e., 
0.1% of the mass of the galaxy). The latter value is consistent 
with 0.14% ± 0.04% the mass fraction in a central S MBH that 
is expecte d to be found in most elliptical galaxies (Haring & 
|Rix|2004) >. The potential and forces due to the central black 
hole are those of a softened point mass with softening length 
e=10- 5 a. 

The orbital structure of the stationary versions of triaxial 
Dehnen models with various cusp slopes and a range of black 



hole mass fractions have been previously studied (M erritt & 
Fridman][T996} |Valluri & Memttl[T998| |Siopis & Kandrup 
2000). "In this paper we restrict ourselves to presenting mod- 
els with 7 = 1 for the following reasons. First, the most likely 
candidates for triaxial elliptical galaxies with slow figure rota- 
tion are th e more luminous "slow rotators" with shallow cen- 
tral cusps (Emsellem et al. 2007). Second, a Dehnen model 
with 7=1 is quite similar (at least in the inner regions) to 
the density profiles of cosmological dark matter halos with 
the main difference being that cosmological density profiles 
fall off more slowly at large radii (r~ 3 compared to r -4 for the 
Dehnen model). Models with 7 = 2 were studied (but are not 
presented) since their dependence on pattern speed is qualita- 
tively similar to that of the 7=1 models with an SMBH. (We 
note that while our choice of triaxial potential is representative 
of a triaxial elliptical galaxy, it is an incomplete representation 
of a dark matter halo which could have a significant potential 
contribution from a stellar disk.) 

At the present time only a few measurements of the pat- 
tern speeds of figure rotation in triaxial dark matter halos and 
early-type galaxies exist. The pattern speeds of fast ba rs mea- 
sured by applying the Tremaine- We inberg method ( jTremaine] 
|& Weinberg 1984; Meidt et aQ2008| l typically indicate that 
the ratio of the corotation radius to the length of the bar 
Rn/a = [1,1 -4] (for a bar of semi-major axis length a) (e.g., 
Debattista et al. 2002 ||A.guerri et al.|2003l|Corsini]2010t|Bu> 
ney & Tremaine| |2008| l. Hereafter we shall use the quantity 
Rn/cLi to compare the pattern speeds of orbits launched from 
different radial shells (of semimajor axis length a,-) to each 
other and to the pattern speeds of bars. 

The pattern speed of early-type galaxy NGC 2974 has been 
measured b y fitting the pro perties of three rings of recent star 
formation (Jeon g et al.|2007| . This galaxy is normally clas- 

4 Following current convention (e.g., BT08) we use the right-handed screw 
rule (with positive angular momentum vector pointed up) for figure rota- 
tion and a right-handed coordinate system. Note that some previous authors 
jSchwar zschild 1982) used a left-handed screw rule for figure rotation and a 
right-handed coordinate system. 
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sified as E4, but the authors argue that all three rings can 
be accounted for as occurring at resonances if this galaxy is 
a SO galaxy with an extended stellar bar and pattern speed, 
ftp ~ 78 ± 6 km s" 1 kpc" 1 . This pattern speed is only slightly 
slower than that of fast bars. Figure rotation of a triaxial dark 
matter halo was proposed as the cause of the extended spiral 
arms in the blue compact dwarf ga laxy NGC 2915 ( [Bureau 
et al. 1999 Bekki & Freeman 2002 ) whose modeling required 
a halo pattern speed of £l p = 7 ± 1 km s" 1 kpc" 1 . This pattern 
speed is an order of magnitude larger than the maximum pat- 
tern speed measured for dark matter halos produced in cos- 
mological iV-body simulations (il p = 1.01/z km s" 1 kpc" 1 ) and 
nearly two orders of magnitude larg er than the median pattern 



speed (ftp = 0.148/i km s 1 kpc ') (Bailin & Steinmetz 



2004). 



Thus, at the present time, observational constraints on the pat- 
tern speeds of triaxial halos and elliptical galaxies are quite 
uncertain. Our fastest pattern speed tt p = 34 km s" 1 kpc" 1 (for 
Rn/a = 2) is slightly slower than the speeds observed in fast 
bars, our slowest Vt p = 0.58 km s" 1 kpc" 1 (for Rn/a = 40) is 
comparable to fastest measured in simulations of dark mat- 
ter halos. We explore a range of values for £l p between these 
limits. 

The triaxial Dehnen models studied here rotate about the 
short axis (z). Simulations of collisionless dark matter halos 
show that the angular momentum a xes have a mean misalign - 
ment of ~ 25° with the minor axis (Bail in & Steinmetz 2005 ). 
However, the dissipative collapse accompanying galaxy for- 
mation is likely to induce angular momentum transport and 
a higher degree of alignment between the spin axis and the 
short axis of the galaxy. The pattern speed of the triaxial fig- 
ure is given in terms of the "corotation radius", hereafter /Jo- 
in a nearly axisymmetric potential, the corotation radius is the 
radius at which the frequency O c of a closed (almost circular) 
orbit in the equatorial (x-y) plane of the potential is the same 
as the pattern frequency (generally called "pattern speed") Sl p : 



R=R n , z =0 



(9) 



Working in the equatorial plane of our triaxial model we 
set/? 2 = (x/a) 2 + (y/bf. In our models Rq is given in units 
of a, the scale length of the semimajor axis, and ranges from 
Rq = 60 (slowly rotating) to Rq = 2 (very rapidly rotating) (see 
Table fTh. In this nomenclature the stationary (non-rotating) 
modelhas its corotation radius at infinity and is labeled with 
Rq = oo. We do not discuss the orbital structure of more pro- 
late triaxial structures similar to fast bars with higher pattern 
frequencies (Rq < 2) since they have been previously s tudied 
( [Skokos et a l. 2002a b{ [Harsoula & Kalapotharakos[2009] ). To 
give the reader a better physical appreciation of the pattern 
frequencies implied by the co-rotation radii, we convert our 
model units to physical units in an elliptical galaxy. For a 
semi-major axis scale radius a = 5 kpc and a galaxy mass of 
M = 5 x 10 n M Q , the unit of time for the model is given by 



a 3 
GM 



1.49 x 10" 



a 


3/2 


M 


kpc 




10 U M Q 



-,-1/2 



yr. (10) 



For these parameters Table [Tjgives the pattern frequency (il p ) 
and rotation period T p = 2ir/Vl p (in years) for each value of 
corotation radius Rq that was studied. 

2.2. Orbit start spaces 



TABLE 1 

Pattern frequencies of figure rotation of triaxial Dehnen 

MODELS WITH 7=1, M bh = 



Rn 


n„ 


n; 


T* 
V 


(a) 


(program units) 


(km s kpc 1 ) 


(yrs) 


2 


2.57X10" 1 


33.8 


1.82 X 10 8 


5 


7.45 xl(r 2 


9.80 


6.28 X 10 s 


10 


2.87xl(r 2 


3.78 


1.63 x 10 9 


20 


1.09xl0~ 2 


1.44 


4.29 x 10 9 


40 


3.85 xlO" 3 


0.51 


1.21 x 10 10 


60 


2.12xl0" 3 


0.28 


2.21 x 10 10 


oo 


0.0 


0.0 


oo 



* Fora = 5kpc,M = 5x 10 H M Q 

Following Merritt & Fridman ( 1996 ), the model's mass dis- 
tribution is stratified into 20 ellipsoidal shells dividing the 
model into 21 sections of equal mass. In this paper, we fo- 
cus mostly on orbits with energy equal to the potential en- 
ergy at the point the 8th shell intersects the major axis of 
the model (x = 1.6a). (Note that shell 10 corresponds to 
the half-mass radius of the model). In the rotating poten- 
tial, orbits were launched from the equi-effective-potential 
surface which is analogous to launching all orbits in a non- 
rotating model from an equipotential surface. Thus all or- 
bits launched from a given shell have the same Jacobi integral 
(Ej = E - h \n p x r| 2 ), where E is the total energy of an orbit. 
We study orbits on 6 energy shells for one pattern frequency 
but defer a full investigation of orbits in self-consistent poten- 
tials to a future paper ( Vallu np011| ). 

Initial conditions for the orbits were selected in two dif- 
ferent ways such that orbits from all four major families (the 
box orbits and th ree families of tube orbits) were represented 
(de Zeeuw [T985 1. Although models are rotating, the orbits 
were launched, integrated, examined, and classified only in 
the frame that is corotating with the figure, allowing for a 
more direct comparison with the stationary models that have 
been studied in the past. 

In a stationary potential, box orbits are characterized by a 
stationary point on the equipotential surface. In a rotating po- 
tential there are no true box orbits (i.e., orbits with stationary 
points in an inertial frame). However, orbits that are charac- 
terized by a stationary point on the effective potential surface 
(surface of constant Jacobi integral) look and behave much 
like box orbits in the frame that is corotating with the figure. 
Accordingly, we launch orbits at zero velocity on a regular 
grid on one octant of the effective potential surface to obtain 
box-like orbits. In a non-rotating triaxial model, a tube orbit 
is characterized by a finite angular momentum (about either 
the long or short axis) which oscillates between two values of 
the same sign. Thus the magnitude of the angular momen- 
tum is not a conserved quantity, but the sign of the angular 
momentum of a tube orbit remains constantr] Consequently, 
tube orbits avoid the center. 

For exploring the phase spac e structure of a sta t ionary tri- 
axial model it i s customary (|Schwa rzschild 1993; |Merritt &| 



Fridman 1 996| [van den Bosch et a l. 2008) to launch orbits 
uniformly from the "x-z start space" i.e., from one quadrant 
of the x-z plane with v x = v, = with v y > determined by the 

5 We refer to tube orbits as "anti-clockwise" or "clockwise" depending on 
whether their time average angular momentum vector (J x for ;t-tubes and J z 
for z-tubes) is positive or negative, respectively (in a right-handed coordinate 
system) 
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energy of the equipotential surface, this is adequate for sta- 
tionary models since the properties of orbits with v y < are 
simply obtained by relying on the symmetries of the model. 
However, in a triaxial model with figure rotation about the 
short axis, orbits that are launched with v v > in the x-z start 
space are different from those with v v < 0. Since orbits of all 
major families intersect th e intermediate y-axi s we use instead 
the "Y-a start space" ( |Schwarzschild 1982) where all four 
families of orbits are launched from the y-axis of the model 
with v v = and at an angle a between the starting velocity 
vector (perpendicular to the y-axis) and the x-y plane. The 
initial conditions for an orbit are given by 



x=Q 

v r = V cos a 



y = Y 
V v = 



z = 

v, = V sina, 



where V is the magnitude of the total velocity of the orbit 
at that position determined from potential energy at the start- 
ing point and the Jacobi integral of all orbits on that surface. 
By allowing 0° < a < 360° we obtain tube orbits with both 
clockwise and anti-clockwise motions. Figure [Tj shows the 
Y-a start space with 225 different initial conditions marked 
by their orbital types. (We plot -30 < a < 330 to illustrate 
that the retrograde z-tubes wrap around the top and bottom 
boundaries.) The maximum value of Y is determined by the 
intersection of the equipotential surface with the y-axis. 

The orbital types were determined visually by plotting three 
different projections in coordinate space (x-y, y — z, and 
x-z) and by examining their time-averaged normalized an- 
gular momenta about the three principal axes (i.e., (J X /\J X \), 
(Jy/\J y \), and (J Z /\J Z \)). For a box orbit (or strongly chaotic 
orbit), all three components of normalized angular momen- 
tum are approximately zero. For x-axis tubes and z-axis 
tubes, the absolute values of the x and z components are equal 
to unity (within 0.005%). Note that a weakly chaotic or- 
bit associated with a tube family will have (J X /\J X \) ~ 1 or 
(J Z /\J Z \) ~ 1 over a finite number of orbital periods. In a three 
dimensional potential, Arnold diffusion would cause such an 
orbit to eventually ergodically fill the entire energy surface 
availa ble to it and would therefore n ot be identifiable as a tube 
orbit (TLichtenberg & Lie bermanj 1 992\ . Hereafter, we refer to 
such orbits as "tube-like". Note that this criterion for classi- 
fying orbits into three major families relies purely on the rela- 
tive magnitudes of the three components of the time-averaged 
angular momenta but does not indicate that these orbits are 
regular. 

Figure[T]shows how orbits of different types are distributed 
on the Y-a start space. Squares indicate box orbits, diamonds 
indicate inner long-axis tubes, triangles indicate outer long- 
axis tubes, plus signs indicate short-axis tubes while an aster- 
isk indicates an orbit id entifi ed as stochastic by its high dif- 
fusion rate (see Section |2~3j ). Note that since the triaxial fig- 
ure rotates in the anti-clockwise direction about the z-axis, we 
refer to anti-clockwise z-tubes as "prograde" and their clock- 
wise counterparts as "retrograde" (even when a specific model 
is stationary). We will see in Figure [3j that stochastic orbits 
are frequently found along the boundaries (or separatrix lay- 
ers) between regions occupied by different orbit families. The 
general location of the different regular orbit families in this 
start space will remain roughly the same as the potential is 
perturbed with non-zero cusp slope 7, the addition of a cen- 
tral point mass, or as the pattern frequency, orbital energy, or 
the shape of the figure are altered. 
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FIG. 1 . — Location of different orbit families on the Y-a start space in a 
stationary triaxial Dehnen model with c/a = 0.5, T = 0.58, 7 = 0, and My, = 0. 
The lines mark approximate boundaries between major families. Squares (□) 
denote box orbits, asterisks (*) denote chaotic orbits, diamonds (o) indicate 
inner long-axis (x) tubes, triangles (A) denote outer long-axis tubes, and plus 
signs (+) indicate z-tubes. 



2.3. Numerical Analysis of Orbital Frequencies 

In fully integrable potentials like the Stackel potentials, all 
orbits are confined to A^-dimensional tori, with the num- 
ber of degrees of freedom (dof), i.e., the number of spatial 
dimensions. Motion on the torus is defined in terms of the 
fundamental frequencies w,-, i = l,...,N, the rates of change 
of the corresponding angle variables. Expressed in terms of 
Cartesian coordinates, the motion is quasi-periodic, e.g., 



X(t): 



OO 

■■J2 A * eil,kt > 



(id 



where the Vk & are linear combinations, with integer coeffi- 
cients, of the three fundamental frequencies: 



Vk = niMi +n 2 M2 + Tl3M3 



(12) 



and the A* are the corresponding amplitudes. The same is true 
for regular orbits in arbitrary potentials, i.e., for orbits that 
respect at least three isolating integrals of the motion (e.g., 



Licht enberg & Lieberman|1992| l. 

When a regular orbit is followed for many (~ 100) dynam- 
ical times, a Fourier transform of the trajectory yields a spec- 
trum with discrete peaks. The locations of the peaks in the 
spectrum correspond to the frequencies in Equation (JTTJ 
and can be used to compute the three fundamental frequencies 
and the integer coefficients (wi J t,»2.A:,W3i ) that correspond to 
e ach peak (e.g.,|Binney & Sperge lfl982| l. 

|Laskar| ( |1990| l developed a fast and accurate numerical 
technique ("Numerical Analysis of Fundamental Frequen- 
cies," hereafter NAFF) to decompose a complex time se- 
ries representation of the phase space trajectory of the form 
x(t) + iv x (t). Our own implementation of Laskar's algorithm 
uses integer programming to obtai n the fundament al frequen- 
cies from the spectrum (Valluri & Merritt||1998| hereafter 
VM98). We integrated orbits using an explicit Runge Kutta 
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integrator (DOP853) of order eight (5,3) due to Dormand & 
Prince wit h stepsize contr ol and dense output by Harier and 
Wanner ( Harier et al.|1993] >. The orbit integration, orbit clas- 
sification, and frequency analysis were all carried out in the 
frame that is corotating with the pattern frequency of the fig- 
ure. 

Even for regular orbits, the character of the motion depends 
critically on whether the three fundamental frequencies are 
independent or whether they satisfy one or more nontrivial 
linear relations of the form 



component is given by 



iw\ + mu)2 + nuji = 



(13) 



where (£,m,n) are integers, not all of which are zero. Gener- 
ally there exists no relation like Equation ( [13) ; the frequencies 
are incommensurable, and the trajectoryfills its torus uni- 
formly and densely in a time-averaged sense. When one or 
more resonance relations are satisfied, however, the trajectory 
is restricted to a phase space region of lower dimensionality 
than N. 

In 2 dof systems, a res onance implies a closed, or peri- 
odic, orbit, e.g., a "boxlet" (|Miralda-Escude & Schwarzschild 



1989). In three dimensions, a single resonance relation like 
Equation ( p~3^> does not imply that an orbit will be closed; 
rather, it restricts the orbit to a space of dimension two (|Mer-| 
ritt & Valluri||1999|l. An orbit satisfying one such relation 



is "thin," confined for all time to a (possibly self-intersecting) 
membrane. In order for an orbit in a 3 dof system to be closed, 
it must satisfy two such independent relations; such orbits are 
likely to be much rarer than thin orbits. In what follows, we 
will use the term "resonant" to refer both to thin orbits (satis- 



fying one relation like Equation 13 1 and to closed orbits (sat- 
isfying two such relations), e.g.7~boxlets. Stable resonances 
of both sorts generate new families of regular orbits whose 
shape mimics that of the parent orbit. Unstable resonant tori 
are typically associated with a breakdown of integrability and 
with chaos. 

Once the fundamental frequencies and their amplitudes are 
obtained, it is possible to obtain two complementary represen- 
tations of phase space at a given energy. A frequency map is 
obtained by plotting ratios of pairs of frequencies (e.g., ui x /uj z 
versus u) y /cj z ) for many thousands of orbits in the potential. 
Such a representation of phase space (e.g., Figure |3ta) and 
Figure|2|a)) is particularly useful for identifying the most im- 
portant orbital resonances. Stable resonances appear as filled 
lines with an increased density of orbits clustered along them. 
This is because these orbits have been trapped by the stable 
resonance. In contrast, unstable resonances appear as "blank" 
o r depopulated l ines and are associated with stochastic orbits. 

|Laskar ([1990) demonstrated that since most orbits in realis- 
tic galactic potentials are only weakly chaotic, they lie close to 
regular orbits in phase space mimicking their regular behavior 
over finite time intervals. Consequently, frequency analysis 
can be used to distinguish between regular and chaotic orbits. 
The frequency drift can give a measure of the degree of chaos 
in an orbit. The frequency drift of an orbit can be determined 
from the change in its fundamental frequencies measured over 
two consecutive time intervals. 

In applying Laskar's formalism, we integrated each or- 
bit for 1007b, where 7b is the period of the long-axis or- 
bit of the same energy in the stationary model. The time 
interval was divided into two equal segments labeled t\ and 
t2 and three fundamental frequencies uj x (t\),ujy(t\), uj z (t\) and 
^tfeXw-yfe), Wzfe) were computed in each time segment. The 
"frequency drift" or "diffusion parameter" in each frequency 



log(A/J = log 
log(A/ v ) = log 
log(A/ z ) = log 



UJ.x(tl)-Ux(t2) 



U}y(tl)-U}y(t 2 ) 

uj z (h)-u z {t 2 ) 



u z (h) 



(14a) 
(14b) 
(14c) 



We define the diffusion parameter log( A/) to be the value as- 
sociated with the largest of the three amplitudes A x ,A y , and 
A z . The larger the value of the diffusion parameter the more 
chaotic the orbit. Our second representation of phase space, 
a diffusion map, is obtained by plotting the initial conditions 
of many thousands of orbits and adding a gray scale (or color 
intensity scale) corresponding to the diffusion parameter. In 
a diffusion map (e.g., Figure [2|b) and Figure l3jb)) the gray 
scale (or intensity of the color) corresponds to log(A/) such 
that regions of initial condition space occupied by regular or- 
bits are white and those occupied by chaotic orbits are dark. 

Several properties of the phase space can be inferred from 
frequency and diffusion maps. We begin with a discussion of 
phase space maps for the stationary (here after "box orbit") 
start space for a baseline model, namely a triaxial Dehnen 
model with c/a = 0.5, T = 0.58,7 = 0, and M bh = 0. The fre- 
quency map in Figure |2ja) shows a plot of the ratios of fun- 
damental frequencies oj x /oj z versus (jJ y /u> z for each of 9408 
box orbits dropped with zero velocity from the equipotential 
surface corresponding to the 8th shell. 

Note that although we focus primarily on the behavior of 
orbits launched from shell 8th, the behavior of an orbit in re- 
sponse to figure rotation depends both on the (radial) energy 
shell from which it is launched as well as the co-rotation ra- 
dius. However, since the density profile has an approximately 
power-law profile outside the inner cusp, the behavior of or- 
bits depends primarily on Rq / a,, where a, is the semi-major 
axis of the rth shell. Consequently, the behavior of orbits 
launched from an outer shell at a moderate pattern speed can 
resemble that of orbits launched at an inner radius and a faster 
pattern speed. Consequently, we will frequently give values 
of Rfi I a t as well as Rq . 

The frequency map shows that most points at the bottom 
left corner of the plot lie on a fairly regular grid which re- 
flects the regular distribution of initial conditions. Away from 
this region the regular, grid-like structure is disrupted by the 
appearance of dark lines corresponding to stable resonances 
and their associated regular orbits. All the orbits lying along 
such lines obey a single resonant condition like that of Equa- 
tion ( 13 1. In the plot, several such resonances are highlighted 
with dashed lines and are labeled with their defining integers 
(£,m, and n). 

Figure[3|a) shows a frequency map for 10,000 orbits initial- 
ized on the Y - a start space with the energy corresponding to 
the x-axial orbit on shell 8 in this model. In this frequency 
map, all four major orbit families are represented: box orbits 
are black, anti-clockwise x-tubes are blue, clockwise x-tubes 
are cyan, clockwise z-tubes are ochre, and anti-clockwise z- 
tubes are red. Note that a significant fraction of the z-tubes 
(red and ochre points) and x-tubes (blue and cyan points) ap- 
pear clustered along straight lines. These are orbits that lie 
close to the thin-shell tube orbit parent of the family and ap- 
pear as "resonance lines" in the frequency map although they 
are not traditionally viewed as resonances. To accurately rep- 
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FIG. 2. — Two representations of phase space for 9408 box orbits launched from shell 8 in a non-rotating triaxial Dehnen model (c/a = 0.5, T = 0.58, 7 = 0, 
and Mti, = 0). (a) A frequency map: d x /u) z versus u) y /u> z for the fundamental frequencies given by the NAFF algorithm. Dashed lines mark various important 
resonances labeled by their integers (i,m, and n). (b) Diffusion map: gray scale corresponds to the diffusion parameter (log(A/)) of orbits at various initial 
positions on one octant of the equipotential surface (the "stationary (box orbit) start space"). Dark regions on the map correspond to chaotic orbits and white 
regions correspond to regular orbits. The labels "X", "Y", and "Z" mark the intersections of the equipotential surface with the three principal axes. 




FIG. 3. — (a) Frequency map and (b) diffusion map for 10,000 orbits in the "Y-a start space" for a model with c/a = 0.5, T = 0.5, and 7 = 0. Orbit families are 
color coded so that box orbits are black, anti-clockwise jc-tubes are blue, clockwise x-tubes are cyan, clockwise (retrograde) z-tubes are ochre, and anti-clockwise 
(prograde) z-tubes are red. In the diffusion map only the chaotic orbits appear colored since the intensity of the color depends on the orbital diffusion parameter. 



resent the resonances in the two tube orbit families it is essen- 
tial to obtain the fundamental frequencies in cylindrical coor- 
dinates about the appropriate symmetry axis. VM98 showed 
that such a representation reveals several unstable tube orbit 
resonances. These can be seen as dark bands in Figure|3jb). 

Fi gure|2|b) is a diffusion map for 9408 box orbits launched 
with zero initial velocity from the equipotential surface (at 
shell 8) of a stationary maximally triaxial Dehnen model with 
7 = 0. The large swath of regular orbits at the bottom left 



of the diffusion plot consists of regular box orbits that origi- 
nate close to the x-axis. In addition, several white bands mark 
regular islands of resonant boxlet families. Each white band 
corresponds to a dark resonance line in the frequency map in 
Figure [2ja). Most white bands are flanked by narrow dark 
regions (occupied by stochastic orbits) which lie along the 
"separatrix" (transition layer) between different orbit families. 
Stochastic orbits (dark regions) are also seen at the intersec- 
tions of stable resonances (white bands), when the parent pe- 
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riodic orbit is unstable. The prominent dark band of chaotic 
orbits that runs along the right edge of the diffusion map cor- 
responds to orbits originating close to the y-z plane. Orbits 
in this region are chaotic due to th e well-known i nstabilities 
of the y-axial and z-axial orbits fHe iligman & Schwarzschild 
1979 Heisler et al. 1982; de Zeeuw & Merritt 1983; Adams 



et al. 



2007| l as well as due to the instability of orbits confined 
in the y-z plane to perturbations perpendicular to this plane 
( Adams et al.|2007) . This region will be seen to expand with 



increasing figure rotation and is referred to hereafter as the 
"y-z instability band". Most of thes e features were described 
in detail in prev ious papers (Valluri & Merritt| [T998 Merritt 
|& Valluri|1999) . 

Figure |3|b) shows a diffusion map for 10,000 orbits 
launched from the Y-a start space. The orbits were classi- 
fied into three different orbit families based on their time- 
averaged normalized angular momentum values. The boxes 
are colored black, anti-clockwise x-tubes are blue, clockwise 
x-tubes are cyan, clockwise (retrograde) z-tubes are ochre, and 
anti-clockwise (prograde) z-tubes are red. Note that regions 
of more intense color (occupied by stochastic orbits) gener- 
ally appear at the separatrix between the major orbit families. 
Within the regions occupied by x-axis tubes (blue and cyan re- 
gions), the transition between the inner and outer x-axis tubes 
is also marked by a weakly stochastic layer. Thick stochas- 
tic layers separate the z-tubes from the box orbits that lie at 
values of Y < 0.5. In addition, stochastic orbits are found at 
a ~ 90° ± 30° and a ~ 270° ± 30° mostly at very small Y val- 
ues (these orbits are launched nearly along the z axis which is 
known to be unstable). 

The main features identified on the frequency and diffusion 
maps above, for the 7 = and My, = case, are meant to give 
readers insight into how they should interpret the results for 
models with figure rotation. 

3. RESULTS 

We now discuss the effects of increasing the pattern fre- 
quency of figure rotation on box-like and tube-like orbits us- 
ing the phase-space mapping techniques discussed in the pre- 
vious section. 

3.1. Effects of slow/Moderate pattern frequency on boxlike 

orbits 

Box orbits were integrated, from start spaces such as the 
one in Figure |2|b), for various model parameters and in- 
creasing pattern speeds. We generally present results for 
a default model which is close to maximally triaxial, with 
c/a = 0.5, T = 0.58 and central density cusp of 7 = 1. We also 
present results for models with a central point mass represent- 
ing a central SMBH with M bh = 0.001 (0.1% of the galaxy 
mass). (Models with 7 = 2 were studied but are not pre- 
sented since their dependence on Rq is qualitatively similar 
although they have a higher fraction of chaotic orbits, VM98. 
Higher black hole mass fractions were previously studied by 
VM98 and were found to induce more chaos.) In each fig- 
ure, we show phase space maps for a non-rotating model 
(Rfi = 00) and three additional values of the corotation radius 
(Rn = 40, 10, and 5). In all cases orbits were launched with 
energy E equal to that of the x-axial orbit started from the 
8th shell. Therefore, for an orbit launched from the 8th shell, 
these pattern speeds correspond to Ru/a, = 24.8,6.2, and 3.1. 

Figure HI shows phase space diffusion maps for box orbits 
in our default model with no central point mass. The non- 
rotating model (panel (a)) is similar to (Figure |2jb)) but con- 



tains a somewhat larger fraction of stochastic orbits. Almost 
all the regular orbits (white regions on the map) are now as- 
sociated with resonant islands containing boxlet families. All 
orbits that qualitatively resemble box orbits in cored poten- 
tials are mildly to strongly ch aotic. 

As mentioned in Section[T| Gerhard & Binney (1985) first 
proposed that the box orbits rendered chaotic due to scattering 
by a massive central point mass could be stabilized in a rotat- 
ing frame due to the envelope doubling effect of the Coriolis 
force. Since a box orbit (in a stationary model) has no net 
angular momentum it oscillates between stationary points re- 
versing its sense of progression around the center each time 
it reaches a turning point. In a rotating frame this means that 
the path described during the prograde segment of the orbit is 
not retraced during the retrograde segm ent because the Cor io- 
lis force on the two seg ments differs ( |Schwarzschild| 1 982] [de] 
|Zeeuw&Merritt|1983| l. 

We see in Figure 141 that this prediction does not hold up 
in realistic triaxial galaxy models. In fact, as the pattern fre- 
quencies of figure rotation increase (Rq decreases) the area of 
the diffusion map occupied by regular orbits (white) decreases 
and the resonant islands shrink until only a small fraction of 
orbits remain regular at R^ = 5. When a central point mass 
is added, the increase in the fraction of chaotic orbits is even 
more significant (Figure [5} and in addition the delineation 
between the various resonant box orbit families is even less 
clear - pointing to increased resonance overlap. We see from 
a comparison of the Figures 4 and 5 that high pattern speed 
increases the fraction of chaotic orbits more significantly than 
the presence of a central black hole. 

We note in passing that a new regular orbit family appears 
as a white region in the lower right hand corner of the dif- 
fusion maps in Figures |4] and |5] for Rq = 10,5. In the sta- 
tionary model, box orbits originating here are highly unsta- 
ble since they lie close to the unstable y-axis. Figure [TTJ top 
row) shows two projections of one such orbit in a modeTwith 
Rn =5. Although it is launched with no net angular momen- 
tum in the rotating frame, this orbit is a short-axis tube. This 
tube orbit (like others in triaxial potentials) does not conserve 
angular momentum, but its value oscillates between two val- 
ues (in this case zero and a large negative value) indicating 
that its motion in the rotating frame is retrograde to that of the 
figure. This new family of orbits lies close to the x-y plane. 
This family appears to arise when the pattern frequency of 
the figure is high enough that Coriolis and centrifugal forces 
cause it to loop around the center. Orbits in this region of start 
space were chaotic in the stationary mode l and are indeed sta- 
bilized by figure rotation as predicted by Gerhard & Binney 
(1985). However, this stabilization results from a complete 
transformation of the orbital character from box-like to tube- 
like, rather than due to a small deflection of the chaotic box 
orbi t aro und the destabilizing center. We will see later (Sec- 
tion |3.3| l that as the corotation radius moves inward, this re- 
gion of stable tube-like orbits grows steadily till most of the 
box start space is occupied by regular orbits, similar to this 
one. 

FiguresHJa) and[5ja) show that most (but not all) the regular 
orbits in the stationary model, which appear as white regions 
on the diffusion maps, are associated with resonant orbits. 
There are two reasons why the majority of box orbits (which 
are associated with two-dimensional resonance surfaces) are 
not stabiliz ed by slow to moderate figure rotation as predicted 
b yJGerhard & Binney| ( [T983] l. 

First, the increase in the fraction of chaotic orbits can be 
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FIG. 4. — Diffusion maps of orbits launched from the box orbit start space at 8th energy shell in models with 7 = 1 and M/,/, = and varying pattern frequencies 
as indicated by labels. The value of Rq = oo refers to a model with no rotation, while Rq = 40, 10, and 5 refer to the corotation radii of models with increasing 
pattern frequencies. These pattern speeds correspond to So/a, = 24.8,6.2, and 3.1. 



attributed to the distinct configurat ion space structure of res- 
onant orbits. As described by |Merritt & Vallu"n] ( 1999) , the 
parent of a resonant orbit family satisfies a condition like 
Iuj x + mWy + nu) z = 0. This additional resonant condition con- 
fines a resonant orbit to the surface of a two-dimensional torus 
in three-dimensional phase space. In general, such orbits are 
not periodic. The more commonly known (but less numer- 
ous) ne ar-periodic versions of such resonances w ere identi- 
fied by Miralda-Escude & Schwarzschild (1989) and given 
names such as "bananas", "pretzels", "fish etc. based on 
their projected shapes. A non-periodic resonant orbit there- 
fore, in comparison, appears similar to a three-dimensional 
box orbit in projection but occupies only a two-dimensional 
surface in configuration space. Consequently, all families of 
stable resonant orbits are centrophobic and avoid the center of 
the potential due to their thin sheet-like structure. The parent 
of the resonant family is surrounded by a region occupied by 
nearly resonant orbits which are characterized by two of the 
same frequencies as the resonant parent, but with an increas- 
ing third frequency. As the third frequency increases from 
zero, the thickness of the orbit increases and the nearest dis- 
tance of approach to the center of the potential (the pericenter 
distance) decreases. When the orbit becomes thick enough to 
pass through the center it is destabilized. The instability at the 



center is a consequence of the divergent central density cusp 
or central black hole present in all realistic galaxy models. 

In the rotating frame of a slowly rotating model, the Cori- 
olis force can produce "envelope doubling" which broadens 
a nearly resonant box orbit driving it into the divergent cen- 
tral cusp (or central black hole), reducing the sizes of stable 
resonant islands, and eventually destroying them. The effect 
of envelope doubling can be seen in Figure [6] which shows a 
resonant boxlet orbit associated with the (3,-1,-1) resonance. 
The top row shows two Cartesian projections of the orbit, x-y 
(left) and x-z (right). In projection (top row) the orbit appears 
similar to three-dimensional box orbits. However, the sec- 
ond row shows intersections of the orbit with the x-y plane 
(left) and the x-z plane (right) shows its perfectly sheet-like 
structure with a clear "hole" in the center. The next three 
rows show cross sections of the same orbit in models with 
Rn = 40,20, and 10. When R n = 20 (4th row) the cross sec- 
tion shows that the Coriolis force broadening has caused it 
to nearly fill the center in the x-z projection. When Rq = 10 
(5th row) the orbit is thick enough to pass through the center 
rendering it completely chaotic. 

While the majority of the regular orbits on the stationary 
start space in Figure 14] and Figure [5] are associated with res- 
onances, a few regular box orbits are non-resonant. These 



10 



Deibel, Valluri, & Merritt 




FIG. 5. — Diffusion maps of orbits launched from the box orbit start space on the 8th energy shell of a model with cja = 0.5, T = 0.58,7 = 1, and M/,1, = 0.001. 
The value of = oo refers to a model with no rotation while Rq = 40, 10, and 5 refer to the corotation radii of models with increasing pattern frequencies. 



too appear to be destabilized by figure rotation. Additionally, 



(Gerhard & Binney 1985 ) predicted that orbits that are chaotic 
in the stationary frame should become regular in the rotating 
frame, but this was not seen. 

The second reason for the increase in the fraction of chaotic 
boxlike orbits can be seen by plotting frequency maps. Fre- 
quency maps for the default model (c/a = 0.5, T = 0.58,7 = 1) 
and Mbh = 0) show resonant orbits clustered along resonance 
lines, and non-resonant regular orbits and stochastic orbits 
scattered over the rest of the map (Figure |7|a)). The points 
scattered along a line with a slope of approximately unity on 
the right hand side of the map correspond to chaotic box orbits 
which lie at the separatrix boundary with a family of short- 
axis tubes with oj x ~ ui y . 

For the slowest rotation frequency (Rq = 40), we observe 
the disappearance of some of the weakest resonances but the 
overall structure of the frequency map remains intact. As the 
corotation radius of the model decreases (higher rotation fre- 
quencies) the entire frequency map appears to shrink with 
the boundaries moving toward the bottom right toward the 
short-axis tube "resonance". This is because as the pattern 
frequency of the model increases, each orbit experiences a 
centrifugal force that changes both uj x and ui y . For a small 
fraction of orbits that lie close to the equatorial plane of the 



model with large y values, the change in orbital frequencies 
is large enough to convert them to resonant z-tubes (e.g., Fig- 
ure 1 1 (top row)). 

The frequency maps show the increased width of the stable 
resonance lines (which occur due to the envelope doubling 
discussed above). In addition a few new resonances appear as 
the pattern speed increases. The shrinking of the frequency 
map, the broadening of resonances, and the appearance of 
new resonances together contribute to a significant increase in 
the overlap of resonances. Resonance overlap is a well known 
cause of global chaos in Hamiltonian systems (Chirikov 1979) 
and may be thought of as occurring when several different res- 
onances compete to trap the same orbit (BT08). This implies 
that even if figure rotation could a llow chaotic box o rbits to 
avoid the center as suggested by Gerhard & Binney (1985), 
the compression of the range of frequencies resulting from 
the modulation of orbit frequencies with the pattern frequency 
results in such significant overlap in resonant orbits that the 
effects of global chaos dominate the behavior of box orbits. 

3.2. Effect of moderate figure rotation on tubelike orbits 

We now explore the effect of figure rotation on tube-like 
orbits in the default model. The diffusion maps in Figure [8] 
show orbits launched from the Y-a start space, with the gray 
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FIG. 6. — Resonant orbit associated with the (3,-1,-1) resonance. The top 
row shows two Cartesian projections of the orbit. The next four rows show 
cross sections of the orbit with the x-y plane (left) and with the x-z plane 
(right) with increasing pattern frequencies (Sq = oo,40,2(), and 10 which 
correspond to Rn/cij = 24.8, 12.4, and 6.2). 



scale representing the diffusion parameters. (Compare with 
Figure[3]for the locations of the major orbit families.) 

As the figure rotation of the model increases, the bi- 
symmetry about a ~ 150° is broken. For Rq = 5, the pro- 
grade z-tubes shrink to occupy a region a ~ 180° ±20° (i.e., 
close to the equatorial plane of the model) but now extend to 
smaller Y values than in the stationary model (panel (a)). In 
contrast, the retrograde z-tubes now occupy a large range of 
a values, but predominantly at large Y). As the co-rotation 
radius decreases, the regions occupied by long-axis (x) -tubes 
also shrink. A greater portion of phase space at small values 
of Y associated with box orbits becomes chaotic (consistent 
w ith what was se en in Figure [4). 

Binney ( 1981) found from stability of the analysis of closed 
periodic (1:1) orbits in the equatorial plane that as figure ro- 
tation increases, retrograde closed orbits become unstable in 
an annular region (the "Binney instability strip") where the 
perpendicular frequency is reso nant with the epicycl ic fre- 
quency. Based on this finding JSchwarzschiId| ( |1982| l erro- 



neously concluded that retrograde tube orbits would become 
increasingly unstable with increasing figure rotation and con- 
sequently chose to largely omit them from the self-consistent 
models he constructed. 

The decrease in the overall fraction of initial condition 
space occupied by the prograde z-tubes and the correspond- 
ing growth of the retrograde z-tubes appears to be in conflict 
with this expectation. However, we do find that at the high- 
est pattern frequency (panel (d)), prograde z-tubes close to the 
equatorial plane (a ~ 180°) are stable over a wider range of 
Y values than the retrograde variety (a ~ 0°). 

In the stationary model (panel (a)) the x-tubes are found 
symmetrically about launch angles of a ~ 90° ± 30° , 270° ± 
30° (orbits are launched perpendicular to the x-y plane). How- 
ever, as the pattern frequency of figure rotation increases the 
portion of start space occupied by x-tubes decreases and or- 
bits with a < 90 and a > 270 are replaced by retrograde z- 
tubes (at Y > 0.6) and box (and chaotic) orbits (at Y < 0.6). 
Thus, in the more rapidly rotating models any orbit launched 
vertically (with no angular momentum about the z-axis) or 
with velocity opposite to the direction of figure rotation will 
"fall behind" the model. The orbit can become a retrograde 
z-axis tube (in the frame corotating with the figure), a box or- 
bit, or a chaotic orbit. (Note that although some orbits are 
"retrograde" in the corotating frame, they could be rotating 
prograde when viewed from an inertial frame.) 

Note that as the corotation radius decreases, inner long-axis 
(x) tubes launched close to the center of the model (small Y 
values) become strongly chaotic, while those launched with 
intermediate Y values are transformed into outer long-axis 
tubes and exhi bit th e "anomalous" behavior discussed by 
1982) (i.e., they are tipped about the y-axis 
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due to the Coriolis forces). The destabilization of inner long- 
axis tubes is also a consequence of the envelope doubling that 
arises due to the Coriolis forces. Figure [9] shows projections 
as well as planar cross sections of two inner long-axis tubes. 
The two columns on the left correspond to an orbit with a 
small pericenter radius which is destabilized, while the orbit 
on the right goes through a phase where it appears chaotic but 
is once again stable at a higher pattern frequency, although it 
is tipped about the y-ax is and significantly flattened. 

Valluri et al. ( 2010| l showed that in triaxial A^-body halos 
x-tubes are the second most important family after the boxes 
(~ 85%) and constitute 11% of the orbits. They also showed 
that in prolate dark matter halos (the shapes most com monly 
found in cosmological simulations; |Dubinski & Carlberg] 
[1991] |Jmg & Suto|[2000l |Allgood et 'al ||2006| > the fraction 
of long-axis tubes can be as high as 75%. Our finding that 
x-tubes with small pericenter radii are easily dest abilized is 
consistent with observations by Valluri et al. ( 2010[ ) that long- 
axis tubes with small pericenter radii are quite easily destabi- 
lized by the presence of a central point mass in self-consistent 
(but non-rotating) A^-body models. 

The discussions in this section indicate that as figure ro- 
tation increases the fraction of phase space occupied by pro- 
grade z-tubes and x-tubes (both clockwise and anti-clockwise) 
decreases and this decrease is accompanied by an increase in 
the fraction of stable retrograde z-tubes and chaotic orbits. 

3.3. Effects of Fast Rotation 

We briefly explore the behavior of orbits when the pattern 
frequency of figure rotation approaches that of bars. Most 
bars are rapidly rotating with corotation radii roughly 1—1.4 
times the length of the bar (e.g., Debattista et al. 2002 ; Aguerri 
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FIG. 7. — Frequency maps for box orbit start space at shell 8 in the default model (c/a = 0.5, T = 0.58, 7 = 1, and M^, = 0) for a non-rotating model = oo) 
and three rotating models (Rn = 40, 10, and 5). 



|et aL|2003|rCorsirii|2010|JBiDney & Tremaine|2008) . 

Figure [TO] shows the effect on orbits launched from the 8th 
energy shell in the default model with figure rotation corre- 
sponding to a corotation radius Rq = 2 (Rn/a, = 1.24). Panel 
(a) shows the diffusion map for orbits launched from the sta- 
tionary start space in the default model; in contrast with Fig- 
ure Qd) (Rn = 5), which was occupied almost entirely by 
chaotic boxlike orbits, at a high pattern frequency the diffu- 
sion map shows a transition to predominantly regular orbits. 
An inspection of the corresponding frequency map (panel (b)) 
reveals that the orbits no longer have box-like characteristics 
but instead predominantly lie along a diagonal line indicating 
that they are more like z-tube orbits. Panel (c) of Figure [T0| 
shows that while a significant fraction of orbits are chaotic, 
the major tube families (the x-tubes and z-tubes) do persist, 
and stable retrograde z-tubes are again a dominant population. 

Examples of some of the regular orbits launched from the 
box start space that appear at the highest pattern frequency 
are shown in the lower three rows of Figure [TT| While these 
orbits do loop around the z-axis and therefore superficially 
resemble z-tubes they are in fact a distinct family more closely 



resembling orbits in barred galaxies. This is not surprising 
since this pattern speed, Ra/at = 1.24, is in the range of values 
seen in fast bars. 

As mentioned previously, the normalized angular momen- 
tum (J Z /\J Z \) is close to unity for bona fide z-tubes. However, 
these new loop-like orbits have a wide range of (J Z /\J Z \) val- 
ues, most of which are significantly smaller than unity. This 
is because although these orbits loop about the center on aver- 
age in a retrograde sense, these orbits frequently make small 
prograde "epicycle-like" loops (instead of turning points in 
the stationary frame) (see left-hand plots of Figure 1 1 For 
these orbits (J Z /\J Z \) is larger than the other two components 
of angular momentum, but it is smaller than unity. Our auto- 
matic orbit classifier therefore identifies these orbits as boxes 
and not z-tubes. These orbits are not axisymmetric in the x-y 
plane and have elongated box-like projections in the x-z plane 
which would be useful in the construction of self-consistent 
triaxial models. 
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FIG. 8. — Diffusion maps for orbits launched from the Y-a start space in a Dehnen triaxial model (with c/a = 0.5, T = 0.58, 7 = 1, and My, = 0). Panel (a) 
shows the model with no rotation (R^ = oo) while the other three panels show models with increasing pattern frequency (Rq = 40. 10, and 5). 



3.4. Dependence on radial shell 

In order to draw conclusions about the orbit composition 
and stability of equilibrium triaxial galaxies with varying pat- 
tern speeds it is not adequate to restrict the study to orbits with 
a single energy (e.g., the 8th energy shell). It is important to 
assess how the fractions of various orbit families depend on 
their potential energy, in self-consistent models such as those 
constructed via Schwarzschild's orbit superposition method. 
Self-consistent Schwarzsc hild models will be constructed in 
a future study ( Vallu rTpOl l| l. In the next two subsections, we 
will study the properties of orbits at six different radii in a 
stationary model as well as in a rotating model with Rq = 10. 
Such an assessment, while not a substitute for constructing 
self-consistent models, can lead to insights into properties of 
such models. We consider orbits in the default model for the 
"Y-a" start space as well as orbits from the "stationary" start 
space. For the latter we consider only a model with a central 
point mass (representing an SMBH). Orbits launched from 
start spaces equipotential surfaces corresponding to shells 2, 
5, 8, 10, 12, and 16 (recall that outer edge of the 16th shell 
encloses 80% of the total mass of the model). These shells 
have major axes (in units of a) of 0.45, 0.95, 1.61, 2.22, 3.10, 
and 6.87, respectively. For Rq = 10, the corresponding values 
ofRn/ai are 22.22, 10.52,6.21,4.50,3.22, and 1.46. For the 



16th shell, the value of Rn/oi = 1.46 is approaching the value 
of that seen in fast bars, accounting for the increased fraction 
of the new family of regular loop-like orbits. 

3.4.1. Tubelike orbits 

Figure [12] shows diffusion maps of the Y-a start space 
for the six different shells. The maps are color coded as in 
Figure |3jb) with black denoting box-like orbits, cyan, and 
blue denoting clockwise and anti -clockwise .x-tubes, ochre de- 
noting clockwise (retrograde) z-tubes and red denoting anti- 
clockwise (prograde) z- tubes. 

The first five shells show many of the characteristics that 
we saw previously: box orbits and chaotic orbits dominate at 
small Y values in all shells. In shell 2, box orbits and inner 
long-axis tubes are found within the inner half of the diffu- 
sion plot and are often mildly chaotic. In the outer parts of 
the model (shell 16), the box orbits only occupy a small frac- 
tion of the phase space and the tube orbit families (clockwise 
and anti-clockwise varieties) occupy equal areas of the initial 
condition space. As before the various orbit families are sep- 
arated by chaotic separatrix layers and in addition unstable 
resonances appear as dark bands in the regions occupied by z- 
and jc-tubes. 

Figure [13] shows the effect of figure rotation on the orbits 
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FIG. 9. — Example of two inner long-axis (x) tube orbits. The top row shows two projections of the orbit in the x-y (left) and x-z (right) planes in the stationary 
model. The second row shows cross sections of the same orbits with the x-y and x-z planes. The next three rows show that the cross sections become broadened 
as the pattern frequency increases, becoming chaotic and/or tipped about the y-axis for higher pattern frequencies. 



in Figure [12] Note that the corotation radius of Rq = 10 was 
selected because this corotation radius lies outside shell 16 
ensuring that none of the shells experience the fast pattern 
speeds seen in Section 3.3 Only ~10% of the total mass of 
the model lies outside a corotation radius of = 10. 

Figure rotation produces roughly the same qualitative ef- 
fects on orbits in all six radial shells although the details de- 
pend on radius: (1) the fraction of chaotic (and boxlike) or- 
bits (black/intense color) increases relative to the non-rotating 
model, especially in the outer shells, (2) the fraction of x- 
axis tubes (blue/cyan) decreases, (3) in shell 8 and beyond the 
region occupied by prograde z-tubes (red) shrinks (there are 
no equatorial plane prograde z-tubes by shell 16 but two new 
bands appear around a - 90°, 270° f\ and (4) the retrograde 
z-tubes (white regions surrounded by ochre bands) dominate. 

To better quantify the effects of figure rotation on the orbit 
populations we computed fractions of orbits in each of the ma- 

6 The new bands of prograde z-tubes that appear at a ~ 100 and a ~ 260 
at the outer most energy shells are classified by our automatic classifier as 
z-tubes but are in fac t x-tubes that have b een tipped about the y-axis by nearly 
90° as predicted by IHeisleret al.|1982l. 



jor families. As discussed previously, we use time-averaged 
normalized angular momentum values to distinguish between 
boxes, x-tubes, and z-tubes. In addition, we visually inspected 
225 uniformly distributed orbits in each shell to confirm these 
automatic classifications and to also allow us to distinguish 
between inner and outer x-tubes. The results of the orbit clas- 
sification are plotted in Figure[l4]for the models in Figure [12] 
(Rq = oo) and Figure 13 (Rq = 10). In the absence of figure 
rotation (left) the fractions of prograde and retrograde z-tubes 
are exactly the same so the (x) and (+) symbols overlap; at 
a higher pattern speed (right), the ochre curve (retrograde z- 
tubes) and red curve (prograde z-tubes) separate clearly show- 
ing the dominance of the retrograde family, especially at large 
radii. In this plot we do not distinguish between chaotic and 
boxlike orbits, since to do so would require setting an arbi- 
trary "cutoff" value for the diffusion parameter at which orbits 
would be classified as chaotic. Most of the orbits classified as 
boxlike are in fact chaotic, but a large number of the tube like 
orbits (especially the x-axis tubes) are quite strongly chaotic 
as can be seen from the intensity of the blue and cyan regions 
in Figures[T2]and[T3 
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FIG. 10. — (a) Stationary start space diffusion map, (b) stationary start space 
frequency map, and (c) Y-a start space diffusion map for orbits in the 8th 
energy shell in models with c/a = 0.5, T = 0.58,7 = 1. M bh = 0, and Rq = 2 
(Rn/a, = 1.24). 



FIG. 1 1 . — Box orbits in models with high pattern speeds are transformed 
to regular retrograde orbits that circulate about the z-axis and loop around 
the center. Top row: a retrograde short-axis tube orbit that originates close 
to the y-axis on the equatorial plane that first appears at moderate pattern 
speeds (Rq = 10,5). Rows 2-4: examples of regular "bar-like" loop-orbits 
that dominate the regular orbit population of orbits launched from the "box" 
start space at Rn = 2 (Rn/a 2 = 1.24). 

The figure shows that the fraction of box-like (chaotic) or- 
bits increases slightly in the inner shells but increases signifi- 
cantly in the outer shells (but is as high as 40% in many shells 
in the rotating models). Note that in contrast, the fraction 
of box orbits drops with radius in the stationary model (left). 
The fraction of inner x-axis tubes is not very significant even 
in the stationary model but decreases significantly, disappear- 
ing entirely at larger radii in the rotating model. (While inner 
long-axis tubes do not play a significant role in maximally tri- 
axial models, they are elongated along the major axis of the 
potential and are vital, in addition to box orbits, for maintain- 
ing the triaxial shape in significantly prolate models.) The 
outer long (x) axis tube fraction decreases slightly at most en- 
ergies, but especially at large radii. The fraction of prograde 
versus retrograde z-tubes remains almost unchanged in the in- 
ner two shells, but there is a significant decrease in fraction 
of prograde z-tubes (red curve (+) signs) in the outer energy 
shells and a corresponding increase in the retrograde fraction 
(ochre curve x). We note that the behavior of orbits in shell 8 
in response to figure rotation is broadly indicative of average 
behavior of orbits in different radial shells. Thus, it is rea- 
sonable to use this shell to gauge the average effect of figure 
rotation on the entire model, but not to predict the detailed 
radial structure of such galaxies. 

3.4.2. Box-like orbits 

Since box orbits are generally regarded as the most impor- 
tant family for supporting triaxiality, it is also necessary to 
study the radial dependence of their stability. The stability of 
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FIG. 12. — Diffusion maps on the Y-a start space for orbits at six different energy levels in a non-rotating model with (c/a = 0.5, T = 0.58,7 = 1, and My, = 0). 



box orbits is strongly affected by chaotic scattering by a cen- 
tral black hole. Since SMBHs are now believed to reside at 
the centers of all galaxies with a significant bulge/ellipsoidal 
components, we now consider behavior of orbits launched 
from the stationary (box) start space in models with SMBHs. 

Figure 15 shows diffusion maps for orbits launched from 
six different radial energy shells (as labeled) in the absence of 
figure rotation. There is a large fraction of chaotic orbits at 
all radii. At small and intermediate radii, there are numerous 
small resonant islands. In the outer two shells, there are only 
three larger islands associated with regular orbits and these 
are mostly associated with non-resonant box orbits (although 
some small resonant islands persist). 

Figure 16 shows the effect of figure rotation (Rq = 10) on 
the box-like orbits in Figure 15 In the inner three shells we 



the regular island s of resonant orbits that we previously dis 
cussed in Section 



3.1 



The disappearance of resonant islands 
is particularly notable in the inner four shells where the non- 
rotating model (Figure [T5]l shows numerous small islands. In 
the outer two shells (which are closer to the corotation radius), 
we see the onset of regularity associated with orbits close to 
the y-axis. These orbits have z-tube-like characteristics which 
was previously demonstrated in Figures 11 and[T0| 

Very near the SMBH, within its gravitationalsphere of in- 
fluence, motion in non-rotating triaxial galaxies is essentially 
regular; orbits take the form of perturbed Keplerian ellipses 
that gradu ally precess due to th e torques from the stellar po- 
tential (M erritt & Vasiliev|2010| l. The innermost radial shell 
of our models contains a stellar mass ~ 0.05M ga i s» 50Mbh and 
so lies well outside of the influence sphere. 



see that figure rotation has the effect of reducing the sizes of 





0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 



FIG. 13. — Diffusion maps on the Y-a start space for orbits at six different energy levels in a model with (c/a = 0.5, T = 0.58,7 = 1; and My, = 0). For Rq = 10, 
Rn/at ■ 22.22,6.21,5.98,3.22, and 1.54 for radial shells 2, 5, 8, 10, 12, and 16, respectively. 



3.5. The effect of varying the shape of the model 

The projected shape of a galaxy on the sky depends both on 
its intrinsic shape (i.e., the value of T and c/a) as well as its 
orientation relative to the line of sight. Although it is possible 
to measure only the projected shape of an individual early- 
type galaxy it is possible to infer the properties of the distri- 
bution of shapes from the assumption that each individual in 
a large sample should be randomly oriented. Using surface 
photometry of the apparent shapes of elliptical galaxies it has 
been shown that their intrinsic shapes cannot be explained if 
they are exclusively a population of randomly oriented oblate 
spheroids, but the distribution of apparent shapes is consis- 
tent with their being a population of triaxial ellip soids (Fasano 



& Vio| 1991 
MerriftPW 



i 



Lamb as et al.|[T992 
IVincent & Rydenl |2005 ) 



Ryden|[T992l [Tremblay & 
Observationally, it 



appears that the lower luminosity ellipticals tend to be flat- 



ter and more consistent with disky oblate systems while the 
higher luminosity systems are more likely moderately oblate 
triaxial. Strongly prolate, perfectly oblate systems, or per- 
fectly spherical systems are rare and the more luminous el- 
liptical galaxies are probably oblate triaxial. The models we 
investigated in previous sections are close to maximal ly tri- 
axial and consequently could be somewhat unrealistic. Mer- 
[ntt]([T997]l explored the range of self-consistent (non-rotating) 
triaxial galaxies of various shapes that could be constructed 
exclusively with regular orbits, in Dehnen models with 7 = 2. 
He showed that the allowed range of shapes for which all or- 
bits were regular tended to be mostly oblate and oblate triaxial 
with a small number of nearly prolate shapes allowed. 

It is out of the scope of this paper to do a full exploration of 
the dependence of orbital structure on galaxy shape. How- 
ever, we studied the effect of figure rotation (Rq = 10) on 
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FIG. 14. — Fractions of orbits in major families as a function of radial energy shell number in the default models, left: without figure rotation and right: with 
figure rotation. Squares (□) connected by black lines denote box-like orbits (which are predominantly chaotic boxes); triangles (A) connected by blue lines 
denote outer .v-axis tube-like orbits; diamonds (o) connected by green lines denote inner .r-axis tubes; crosses ( X ) connected by red lines denote prograde z-tubes, 
and plus signs (+) connected by ochre lines denote retrograde z-tubes. In the absence of figure rotation (left) the fractions of prograde and retrograde z-tubes 
are exactly the same so the (x) and (+) symbols overlap; at a higher pattern speed (right), the ochre curve (retrograde z-tubes) and red curve (prograde z-tubes) 
separate clearly showing the dominance of the retrograde family, especially at large radii. 



orbits launched from shell 8 in models with three additional 
shapes, with 7=1 and no black hole. Figure [T7] shows the 
Y-a diffusion maps for four different models: (1) a nearly 
oblate model, (2) the nearly maximally triaxial model studied 
in previous sections, (3) a moderately triaxial model, and (4) 
a nearly prolate model. 

In the nearly oblate model (panel (a)) the fraction of x-axis 
tubes is insignificantly small, being almost entirely replaced 
by z-tubes. The few surviving x-tubes are outer long-axis 
tubes. Box orbits continue to occupy only a small fraction 
of this start space. This is the only rapidly rotating model 
in which there is not a large difference between the prograde 
and retrograde short-axis tube orbit fractions (retrograde frac- 
tion^. 42, prograde fraction=0.35). The maximally triaxial 
model in panel (b) has a larger fraction of chaotic orbits than 
models of any other shape. The moderately triaxial model 
(panel (c)) has a significant fraction of z-axis tubes (espe- 
cially retrograde ones) and a moderate fraction of long-axis 
(x) tubes. Almost all the long-axis tubes are associated with 
anomalous orbits and the few that show characteristics of in- 
ner long-axis tubes are mildly chaotic. Most of the box orbits 
in the model are quite stable. Finally, the nearly prolate model 
(panel (d)) is dominated by retrograde z-axis tubes while the 
region occupied by prograde z-tubes is now entirely occupied 
by chaotic orbits. Box orbits and x-tubes are found in smaller 
numbers and some are moderately chaotic. 

Admittedly, it is dangerous to draw conclusions about the 
shape distributions of self-consistent galaxy models with fig- 
ure rotation, from looking at the populations of orbits in a 
selected number of start spaces. However, our analysis shows 
that orbit populations evolve slowly and systematically within 
an individual model (with varying a, Y, and radial shell), as 
well as with the shape of the figure and the degree of figure 
rotation. In particular we see that in models with a range of 
shapes, stable short-axis tubes with retrograde motion are the 
dominant population in the start space. At fast pattern speeds, 
box orbits are also replaced by retrograde loop-like orbits that 
circulate about the z-axis. The dominance of retrograde z- 
tubes at moderate to fast pattern speeds suggests that these 



orbits are very likely to be important in self-consistent triaxial 
galaxies. 

4. SUMMARY AND DISCUSSION 

We have investigated the effects of figure rotation on the or- 
bital structure of a family of triaxial galaxy models that pro- 
vide an accurate representation of the luminosity profiles of 
elliptical galaxies. Most of our study is restricted to models 
that are close to maximally triaxial with c/a = 0.5, T = 0.58, 
with shallow cusps of 7 = 1. Rotation of the figure was as- 
sumed to be about the short axis. 

Our major results are summarized below. 

1 . As figure rotation increases, the fraction of initial con- 
dition space occupied by prograde z- and x-tube orbits 
(both clockwise and anti-clockwise) decreases. The 
fraction of stable retrograde z-tubes and chaotic orbits 
correspondingly increases. 

2. Almost all regular box-like orbits in realistic triaxial po- 
tentials with central cusps and/or black holes are asso- 
ciated with resonant orbits. The Coriolis forces in the 
rotating frame of a triaxial potential destabilize these 
near-resonant box orbits as well as inner long-axis tube 
orbits, by thickening them enough to drive them into 
the center. The outer long-axis tubes are also slightly 
destabilized by figure rotation but less than the inner 
long-axis tubes. Box orbits are further destabilized be- 
cause the modulation of the orbital frequencies by the 
rotational frequency results in a shrinking in the overall 
range of orbital frequencies, and therefore in significant 
resonance overlap. Orbits that are chaotic in the sta- 
tionary model are not in general stabilized unless they 
occupy very special regions in phase space. 

3. Over a wide range of pattern speeds there is little evi- 
dence that orbits that are chaotic in the stationary model 
are stabilized in the rotating model due to avoidance 
of the center. At fairly rapid rotation speeds, a small 
region close to the y-axis that is occupied by chaotic 




box orbits in the stationary model is indeed stabilized 
by fig ure rotation as predicted by Gerhard & Binney 
(1985). However, this stabilization is due to a com- 
plete transformation of the orbital characteristics from 
box-like to tube-like, rather than due to deflection of 
the chaotic box orbit around the destabilizing center. 
At even higher pattern speeds (Figure 10 1 Rn/cij < 3 



box-like orbits begin to show a sudden transformation 
to a new family of retrograde loop-like orbits similar to 
those seen in rapidly rotating bars. These orbits loop 
around the z-axis in an overall retrograde sense, but ex- 
ecute small "epicycle"-like prograde motions at points 
in the orbit that correspond to stationary points in the 
non-rotating model. 

4. A surprising result of our study is that as the pattern 
speed Vl p increases there is a decrease in the over- 
all fraction of prograde short-axis tubes which are in- 
creasingly replaced by the retrograde variety. Prograde 
short-axis tubes only persist when launched at angles 



~ ±20° from the equatorial plane of the model. Orbits 
launched at larger angles to the equatorial plane have 
lower angular momenta about the short axis at a given 
energy, and consequently tend to "fall behind" the fig- 
ure and become retrograde as it rotates past them. 

5. There is a fairly strong dependence of the changes in 
the behavior of orbits on the radial energy level from 
which they are launched. (1) The fraction of inner x- 
tubes decreases slightly at all radii, essentially disap- 
pearing at large radii. (2) The fraction of outer x-tubes 
also decreases at all radii. (3) The fraction of retrograde 
z-tubes increases and is accompanied by a decrease in 
the fraction of prograde z- tubes. (4) There is an overall 
increase in the fraction of chaotic/box-like orbits at all 
radii, but the increase is most significant at intermediate 
and large radii. 

6. As the shape of the triaxial figure changes, the orbital 
populations at shell 8 (roughly the half-mass radius) 
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FIG. 16. — Diffusion maps on the box start space for orbits at six different energy levels in a model with (c/a = 0.5, T = 0.58,7 = 1> an£ i Mm = 0.001) with 
co-rotation radius of Rq = 10. For shells 2,5, 8,10,12, and 16, respectively, this pattern speed corresponds to Rsi/ctj : 22.22, 10.52,6.21,4.5,3.22, and 1.46. 



vary. Yet in all cases with figure rotation (with Rq = 10) 
retrograde z-tubes are the most robust and stable popu- 
lation. The fraction of chaotic orbits is lower than in the 
maximally triaxial models. 



It has been known for three decades, since the work of 
|Schwarzschild| ( |1982[ ), that self-consistent triaxial equilibria 
can be constructed if the triaxial figure is rotating slowly. 
Schwarzschild's study was restricted to galaxies with central 
cores; our study of orbits in more realistic triaxial models sup- 
ports the view that galaxies with slow figure rotation have a 
large fraction of stable "triaxial orbit families" such as boxes 
and long-axis tubes which can be used to construct stable tri- 
axial models. We find, however, that as the pattern speed of 
figure rotation increases, these two important families (sta- 
ble resonant boxes and long-axis tubes) become increasingly 
unstable, due to the effects of envelope doubling. When the 
pattern frequency is between Q p ~ 3-30 km s" 1 kpc"'(for an 



elliptical galaxy of scale radius 5 kpc and mass 5 x 10 n M Q ), 
figure rotation produces a sufficiently large amount of global 
chaos to make it unlikely that equilibrium triaxial models with 
such pattern speeds would exist - at least, if one assumes 
that a significant fraction of elongated regular orbits are re- 
quired. At higher pattern frequencies, comparable to those 
of fast bars, the phase space associated with box-like orbits 
becomes stable once again, although the orbits now resemble 
loops found in fast bars. 

Pattern speeds of real triaxial galaxies and dark-matter ha- 
los are poorly constrained. We explored a range of values 
for figure rotation motivated by theoretical measurements of 
pattern speeds o f dark matter halos ar ising from cosmolog- 
ical s imulations (Bailin & Steinmetz 2004; Brya n & Cress| 
|2007| l as well as one claimed observational measurement of 
rotation of a triaxia l dark matter halo with a pattern speed 
in NGC 2915 (Bureau et al. 1999) and an early-type galaxy 
NGC 2794 HJeong et al.||2007| l. The shape assumed for the 
rigid triaxial dark matter halo potential in the model for NGC 
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FIG. 17. — Effect of varying the shape of the galaxy on the orbit population is seen in Y-ct diffusion maps for models with (7=1 and i?n = 10). Panels 
correspond to models with various shapes as indicated by values of c/a and T. The model (a) is nearly oblate, (b) is nearly maximally triaxial, (c) is moderately 
triaxial, and (d) is nearly prolate. The different orbit families are clearly separated by stochastic separatrix layers. 



2915 ( |Bekki & Freeman |2002) (c/a = 0.6, T = 0.56) is close 
to our default model shape (c = 0.5, T = 0.58). Our study sug- 
gests that the pattern speed of ~ 7 ± 1 km s" 1 kpc _1 required 
to produce the extended spiral features in this galaxy is ade- 
quate to significantly destabilize the two major families that 
sustain the high density along the major axis, namely, the box 
orbits (predominantly associated with resonances) as well as 
the inner long-axis tubes. This does not imply that the spiral 
arms in this galaxy are not generated by tidal torques from a 
triaxial halo but does suggest that the halo is unlikely to be as 
strongly triaxial as assumed by these authors. A rapidly rotat- 
ing prolate-tria xial halo (predicted by c osmologi cal jV-body 
simulations, |Dubinski & Carlb erg 1991] |Jing & Suto||2000"l 
Allgo od et al.J|2006| )) is also unlikely, because inner long 



axis tubes, which are important in such models (jValluri et al. 



2010 1, are easily destabilized by large pattern speeds. 

lithe orbits that are required to construct triaxial galax- 
ies, namely, the boxes and inner long-axis tubes, become 
highly chaotic (or disappear entirely) at intermediate pattern 
speeds, they are less elongat ed than necessary to construct 
self-consistent galaxies (Merritt & Vall uri|1996 i. This places 
additional constraints on both the range of shapes as well as 
the range of pattern speeds of real elliptical galaxies. 



While we did not investigate barred galaxies, our study of 
a few cases with rapid rotation suggests that in addition to the 
limited range of shapes that elliptical galaxies can assume, 
stable galaxies can only have two ranges of figure rotation - 
very slow (with rotation periods T p >5x 10 9 yr) or very fast 
(T p < 2 x 10 8 yr). 

Since some of the orbits in models with rapid figure rota- 
tion can be quite different from those in models with slow 
figure rotation, the implications of our studies could also have 
consequences for galaxy and black hole scaling relations. All 
existing stellar dynamical measurements of SMBH masses in 
barred galaxies (e.g., NGC 1023, NGC 3384, NGC 4151) as- 
sume that the dynamics in the central region is unaffected by 
the presence of a large-scale bar, and consequently these sys- 
tems are modeled as axisymmetric. It has been argued re- 
cently that the "black hole fundamental plane" may in fact be 
an artifact of the dep endence of the M^-a s caling relation on 
host morphology (Graham 2008b ; Hu 2008 ), and in particular 
could be a consequence of the fact that barred galaxies appear 
to be offset from the Mbh - & relation for non-barred galaxies. 
In addition to measurement errors in the velocity dispersion 
parame ter a that arise f rom large scale streaming motions in 
the bar (Graham 2008a), neglecting figure rotation could also 
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lead to errors in Mbk- 

There are very limited observational prospects for measur- 
ing the pattern speeds of individual triaxial galaxies and dark 
matter halos. Measurements have only been possible where 
there is evidence of an extended disk-like component in which 
resonance patterns have been excited by the figure rotation 
or when the Tremaine- Weinberg method is applicable. In- 
creasing the number of observational measurements of pattern 
speeds in individual galaxies is likely to be challenging. Our 
study of the behavior of orbits suggests that triaxial elliptical 
structures may have only two ranges of pattern speeds. 

However, it is possible that pattern speeds of elliptical 
galaxies may be correlated in a statistical sense with observ- 
able measures of angular momentum content such as the Xr 
parameter used to sep arate galaxies into slow and fast rotators 
(Emsellem et al.||2007| ). Certainly, a similar correlation has 
been found between the pattern speeds and spin parameter A 



(Peebles 1969 ) of cosmological dark matter halos (Bailin & 
|Steinm etz 2004 ). Further st udies of self-consistent models 
are in progress ( Valluri|201 1\ and will shed light on both the 
ranges of pattern speeds for which self-consistent equilibrium 
triaxial galaxies exist, as well as the observable properties. 
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